Symmetry of standing waves generated by a point defect in epitaxial graphene 
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Using scanning tunneling microscopy (STM) and Fourier Transform STM (FT-STM), we have 
studied a point defect in an epitaxial graphene sample grown on silicon carbide substrate. This 
analysis allows us to extract the quasiparticle energy dispersion, and to give a first experimental 
proof of the validity of Fermi liquid theory in graphene for a wide range of energies from -800 
meV to +800 meV . We also find evidence of a strong threefold anisotropy in the standing waves 
generated by the defect. We discuss possible relations between this anisotropy, the chirality of the 
electrons, and the asymmetry between graphene's two sublattices. All experimental measurements 
are compared and related to theoretical T-matrix calculations. 
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Graphene is a two-dimensional zero-gap semiconduc- 
tor that shows many fascinating physical properties, in 
particular the quasiparticles near the Dirac point have a 
massless spectrum and linear dispersion. These relativis- 
tic particles have a velocity of about one percent of the 
speed of light, and are able to propagate ballistically over 
mesoscale distances [1[. The two inequivalent A and B 
sublattices of the 2D honeycomb carbon structure give 
rise to the fermionic pseudospin and chirality that lie at 
the origin of the forbidden backscattering in the presence 
of extended impurities [2] , and of the anomalous Friedel 
oscillations in the vicinity of localized impurities [3] . The 
possibility to make few-layer graphene samples by the an- 
nealing of a silicon carbide substrate has been previously 
demonstrated in Refs. [4J, |3, [6J. The Dirac fermion be- 
havior has been evidenced in epitaxial graphene by angle- 
resolved photoemission spectroscopy (ARPES) [7|, [8[. 

Here we analyze the Fourier transform of STM con- 
ductance images in bilayer epitaxial graphene. It has 
been shown that such a study should allow one to extract 
the quasiparticle dispersion. This was first observed in 
superconductors |9[ , then by some of the authors in a 
semi- metallic material [10|, lll| , and more recently in bi- 
layer graphene [12|; the latter experiment was however 
limited to a small energy window (±100meV). Various 
other STM observations of the LDOS distortions in the 
presence of impurities in few- layer graphene exist in the 
literature (see for example 13lJl4JJl5 Ul6j | , for hy drogen- 
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induced defects on graphite surfaces, and [lz 
eral types of defects in graphene bilayers). 

In this Letter we present an in-depth theoretical (sin- 
gle impurity T-matrix) and experimental (STM and FT- 
STM) study of the local density of states (LDOS) in the 
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vicinity of an impurity. We have looked for point de- 
fects and we have focused on a single isolated impurity 
that strongly affects the LDOS of the epitaxial graphene 
layer; the distortion induced by this defect turns out to 
have a threefold symmetry. Using FT-STM we test the 
robustness of the quasiparticle dispersion, and the range 
of validity of the Fermi liquid theory in bilayer (or multi- 
layer) graphene for a large energy window (—800 to +800 
meV). Our experiment was performed with a LT-STM 
from Omicron at 77K in a base pressure of 10 -11 mbar. 
The samples were prepared in UHV by the annealing of 
n-doped SiC(OOOl) at 1000K starting from a 3x3 recon- 
struction [5[. 

Figure [I]A shows a topographic image of a large 
graphene terrace taken at -17 meV (probing full states). 
This layer shows an intriguing "star-like" defect with an 
apparent six-fold (C6v) symmetry. This atomic defect 
is accompanied by a strong distortion of the graphene 
lattice. The center appears black, which is a dramatic 
change from the case of the unperturbed lattice. As 
schematized in figure [T]C, a detailed analysis of the real- 
space image, shows that the point defect is directly lo- 
cated over or under a lattice site. The 2D Fourier trans- 
form of the topographic image is shown in[T]E. Six high- 
intensity regions corresponding to the periodicity of the 
hexagonal lattice appear at the centers T of the six Bril- 
louin zones (BZ) neighboring the first BZ. Six high- 
intensity features (open circles) with an anisotropic in- 
tensity are also observed around the corners of the first 
BZ. They can be attributed to the scattering of quasipar- 
ticles between isocontours surrounding the K points (the 
corners of the BZ) as schematized in figure[Tp. This type 
of scattering is called inter-nodal or inter- valley scatter- 
ing. The six high-intensity regions at the corners of the 
first BZ correspond to scattering processes that link two 
consecutive isocontours. Similar features have also been 
observed in 0, Q E3, • 

Using FFT filtering 



we have removed the lattice 
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FIG. 1: A) STM topographic image (10 x lOnra 2 , -17meV 
In A) showing an isolated defect. B) Zoom-in (5 x 5nra 2 ) or 
the defect, subsequent to a FFT filtering removing the atomic 
resolution and all other features with wavevectors outside the 
yellow circle in E. C) 2.5 x 2.5nra 2 zoom in on B with the 
schematic atomic lattice of graphene (A/black and B/yellow 
atoms) superimposed over the standing-wave pattern arounc 
the point defect. D) schematic representation of the 2D Bril- 
louin zone and of the constant energy contours, indicating the 
possible scattering processes. E) The FFT power spectrum o: 
the 2D topographic image in A. F) The real-space LDOS in z 
bilayer graphene calculated for a point defect using a single- 
impurity T-matrix approximation at 300raeV above the Dirac 
point. 



periodicity vectors and all other features with wavevec- 
tors outside the yellow circle in Fig. [Tp, thus taking into 
account only the intra-valley and the inter-consecutive- 
valley scattering processes (K — ► K\ type) depicted in 
figure [Tp. The resulting real-space image is shown in fig- 
ure [TJ3. This operation strongly enhances the anisotropic 
intensity observed also on the bare topographic image. 
The LDOS near the defect shows a clear threefold (C3v) 
symmetry. While, as depicted in Fig. [T]0, close to the 
impurity we observe a fairly homogeneous standing- wave 
ring, with an almost perfect six-fold symmetry, farther 
away from the impurity, the intensity is clearly higher 
along three axes (drawn in green in Figs. [T]G and [if), 
direction. 

Let us now analyze the possible origin of this three- 
fold symmetry. If the defect is localized above an A 
atom (shown in black in Fig. HP), the nearest neighbors 
are three B atoms (C3v symmetry). The second- nearest 
neighbors are six A atoms (C6v symmetry). Hence, the 
three-fold symmetry comes from the lattice environment, 
in particular from the scattering of electrons between the 
A and the B sublattices. This is related to the pseu- 
dospin and chirality of electrons in graphene and will be 
explained in more detail in what follows. 

The theoretical LDOS modulation for a bilayer 
graphene in the presence of a single impurity is repre- 
sented in Fig. [Tp. This was obtained at an energy of 
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FIG. 2: Calculated real-space LDOS modulations in a mono- 
layer graphene for a point defect placed on top of an A atom. 
Figs. [2] A, B, C depict the contributions of selected scattering 
processes to the LDOS (K —> Ki,K — ► K^ and K — ► K3 
respectively) as indicated in Fig.[lp; in[2F, and[2p we depict 
the separate contributions of the A and B sublattices to the 
K -> K x LDOS modulations in[2K. Figs.[2p and E depict the 
A and B sublattice contributions to the full LDOS, evaluated 
when all scattering processes are considered. In H we depict 
the sum of the A and B contributions to the full LDOS, when 
the weight of the K — ► K\ B component is increased three 
times compared to that of the K — ► K\ A component. 



300meV above the Fermi level, using the formalism pre- 
sented in [20] . One can observe quite a few features sim- 
ilar to the ones depicted in Fig. [T]G, including the exis- 
tence of a three-fold symmetry. However the three-fold 
anisotropy is much less pronounced. 

Figure [2] displays the calculated real-space modulations 
in the LDOS close to the point defect, when only specific 
scattering processes have been considered. For simplic- 
ity this calculation has been performed for a monolayer 
graphene, as for bilayer graphene it is less clear how much 
the two bands and the two layers contribute to the mea- 
sured LDOS. The intra- nodal scattering processes are re- 
sponsible for radially-symmetric features with small if 
wavevectors and should therefore not be responsible for 
the presence of the three-fold symmetry. Hence we do 
not focus on these modulations here. 

In Figs.[2]A,[2B andEJG, we plot the LDOS modulations 
due solely to scattering between nearest, second- nearest, 
and third- nearest neighboring nodes, generically denoted 
K — > K\ , K — > i\~2 and K — > K3 respectively, as schema- 
tized in Fig.[lp. First, it is clear that the most important 
effect of the impurity, which is responsible for the R3 pat- 
terns observed experimentally, is the K — > K\ scattering. 
Second, close to the impurity, the patterns are almost six- 
fold symmetric, with a slight three- fold anisotropy. Thus, 
none of these scattering events contributes significantly 
to the threefold symmetry. To investigate the source of 
the anisotropy, we first note that, in Ref.[3] it was shown 
that the K — > K\ high-intensity regions appearing in the 
FT of the LDOS at the corners of the BZ are anisotropic, 
more pronouncedly for monolayer graphene than for bi- 
layer graphene. This is due to the chirality of the quasi- 



particles in graphene. We believe that this asymmetry 
in the observed momentum space features, and thus the 
chirality of quasiparticles, is responsible for the small 
amount of C3v asymmetry observed in the theoretical 
real-space LDOS calculations. The second observation is 
that the anisotropy comes mainly from the B component 
of the LDOS, that describes the double virtual hopping of 
quasiparticles between the A and B sublattices[2l|. This 
is detailed in Figs. [2f and[2]G, where we plot separately 
the A and the B sublattice contributions to the LDOS 
modulations generated by K — > K\ scattering. We can 
see that the contribution from the A atoms is fully six- 
fold symmetric, while the contribution of the B atoms 
is three- fold symmetric. This is also true for the A and 
B components (depicted in [2)3 and [2f] respectively) of 
the full LDOS obtained from summing the contributions 
of all scattering processes. Fig. [2)H[ depicts the sum of 
the A and B contributions to the full LDOS, where in 
order to increase the anisotropy, we have multiplied the 
K — > K\ B contribution by three. This may mimic bi- 
layer graphene where, due to the coupling between the 
two layers, the sublattices A and B may not contribute 
equally to the observed LDOS. A very strong threefold 
symmetry is observed in this weighted superposition, and 
the resulting image corresponds more closely to the pat- 
tern observed experimentally in figure [IJ>). This is con- 
sistent with having a defect at an A site whose dominant 
effect is the scattering of sublattice B electrons between 
two consecutive valleys (K — > K\). 

We are also interested in how such a defect affects the 
relativistic character of the quasiparticles. To find this, 
we measured the quasiparticle energy dispersion using 
the Fourier transform of conductance images at different 
energies. Every image was acquired using a lock- in am- 
plifier and a modulation voltage of ±20 meV which gave 
the energy uncertainty. We have performed a 2D FFT of 
these images, and in Figs.[3]A and[3^ we show the result- 
ing intensity profiles along the direction K — Y — K. For 
a wide range of energies, between energies well below the 
Fermi level (—950 meV) to about — 150meV, the width 
of the peak at the center of the Brillouin zone (measured 
at half maximum) shows a clear dispersion with the bias 
voltage (it decreases when the energy moves towards the 
Fermi level). At energies closer to the Fermi level the 
profile becomes more complex and displays a dip in the 
intensity profile close to the T point. The central peak 
also shows two lateral structures (shoulders). The width 
of the central feature does not appear to disperse with 
energy close to the Dirac point, but it starts dispersing 
again for higher positive energies. In Fig. [3)0 we show 
the STS measurement of the LDOS (dl/dV as a func- 
tion of bias voltage) taken close to the impurity. We 
can see a 150 meV gap-like feature centered around the 
Fermi level at meV . Fig. [H shows a linear dependence 
of the width of the central feature with energy, except 
in an energy range from —200 meV to 100 meV. For 
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FIG. 3: A) and B) Intensity profiles along the K - Y - K 
direction of the 2D FT of the LDOS taken for different bias 
voltages. C) STS measurement taken near the impurity show- 
ing a gap centered around the Fermi level. 



these energies the shape of the central peak is more com- 
plex. The width of the features at the K points follows 
a similar dispersion. Beyond +500 meV, all the features 
seem to follow a different dispersion branch, as indicated 
in Fig. [H We estimate the Dirac point at —100 meV 
below the Fermi level. The spreading in the dispersion of 
the points close to the Fermi level could be attributed to 
the presence of the gap. We compare the experimental 
results with similar theoretical profiles obtained using a 
T-matrix approximation for a single localized impurity 
in bilayer graphene [3] . No gap at the Fermi level was in- 
cluded in the theoretical calculation, while a small gap of 
^100 meV was assumed near the Dirac point, which was 
taken to be close to —250 meV. In Fig. 2£> we plot the 
dispersion of the central and K-points features, and of the 
shoulder, obtained from the theoretical curves. We also 
observe the presence of two different dispersion branches. 
The second dispersion branch and the shoulder arise be- 
cause of the bilayer graphene upper band which opens at 
energies higher than the inter-layer coupling. For mono- 
layer graphene these extra features should not appear. 

In conclusion, our measurements show that the quasi- 
particle approximation and the Fermi liquid theory are 
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FIG. 4: The dispersion for the width of the central-ring fea- 
ture (red circles) and of the K-points feature (black triangles) 
obtained from the experimental data (A), and from the the- 
oretical K — r — K cuts (B). The dispersion of the central- 
feature shoulder on the theoretical curves is indicated by blue 
squares. 



robust over a large range of energies. The complex struc- 
ture of the central feature (the existence of the shoul- 
der), as well as the presence of two dispersion branches 
are consistent with the bilayer (or multilayer) character 
of the graphene sample. While the point defect modifies 
the electronic wave-function in its vicinity, a clear linear 
dispersion is still observed, and the relativistic character 
of the quasiparticle is preserved. Also, the STS mea- 
surements in Fig. [3)0 indicate the presence of a gap cen- 
tered at the Fermi point, and not at the Dirac point in- 
ferred from our FT-STM measurements. This gap is thus 
different from the Dirac-point gap observed by ARPES 
[3, Isj which was attributed to the different doping levels 
of the epitaxial graphene layers. It is however consistent 
with previous STM measurements performed on epitaxial 
graphene [22|, and more recently on exfoliated graphene 
[23|, where a gap observed at the Fermi level was at- 
tributed to a pinning of the tunneling spectrum due to 
the coupling with phonons. 

A marked three-fold symmetry in the vicinity of the 
point defect was also observed. Detailed theoretical stud- 
ies of the different scattering events underline the impor- 
tance of the scattering between consecutive nodal points 
(K — ► Ki). Our analysis also suggests that for bilayer 
graphene, the two sublattices do not contribute equally 
to the measured LDOS. Thus, the amount of asymme- 
try observed in the LDOS near the impurity is consistent 
with the fact that the scattering of the electrons between 
the sites of the sublattice on which the impurity is located 
(A) and the electrons on the B sublattice dominate over 
the scattering of the electrons between sites of the same 
(A) sublattice. 
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